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We study temperature induced metal-insulator transition in doped ferromagnetic semiconductors, 
described by s-d exchange model. The transition is a result of the mobility edge movement, the 
disorder being due to magnetic ions spin density fluctuations. The electrons are described in the tight 
binding approximation. Using ideas and methods of Anderson localization theory we obtain simple 
formulas, which connect the mobility edge with short-range order characteristics of the magnetic 
subsystem - static spin correlators. We discuss the application of the theory to several groups of 
novel colossal-magnetoresistance materials and include the reproduction of the paper jij published 
by us 10 years ago. 

PACS numbers: 71.28.+d, 71.30. +h, 72.15.Rn, 72.20.-i 



I. PRE-INTRODUCTION 



Conducting ferromagnets exhibiting resistivity peak 
and colossal magnetoresistance near the Curie point 
such as europium oxide- and sullide-based compounds 
EuXi_ y , Eui-^RrX where X=0, S, R=La, Gd and 
chromium-chalcogenide spinels A&2X4 where A=Cd, 
Hg, Cu, Fe, X=S, Se were studied from 60's to late 70's 
extensively. Because of difficulties in sample preparation 
and other reasons the interest of RD community in these 
compounds declined in early 80's. The same, even ear- 
lier, happened to manganite perovskites Lai-xA^MnOa 
where A=Ca, Ba, Sr, Pb. Nevertheless, some studies in 
this field continued through the 80's. In 1993, in connec- 
tion with new thin-film technology for manganite per- 
ovskites H the interest in these, called today colossal 
magnetoresistance (CMR), compounds has been revived. 
As a result, approaches and concepts developed and ap- 
plied during 70's and 80's has again been brought into 
the focus of attention. Partial list includes: magnetic po- 
laron versus FM-AFM phase separation, ferromagnetism 
in narrow conduction band, non-pole structure of Green 
function in s-d and Hubbard models. This concerns also 
the concept of localization driven by magnetic fluctua- 
tions 

We published 10 years ago a series of papers on the 
resistivity peak and magnetoresistive effect in magnetic 
semiconductors WpLQ. To the best of our knowledge, 
the present authors were the first to apply the concept 
of mobility edge for disorder determined by spin fluctua- 



tion. The concept of localization driven by magnetic fluc- 
tuations has become actual again, when the problem of 
resistivity peak and colossal magnetoresistance has again 
been brought into a focus of attention. 

Below we reproduce our paper " Anderson localization 
in ferromagnetic semiconductors due to spin disorder. I. 
Narrow conduction band" jy, the first one from this se- 
ries. 



II. INTRODUCTION 

Magnetic semiconductors form a wide class of materi- 
als with 3 unique physical properties, strong temperature 
dependence of the electroresistivity being one of them [^) . 
This phenomenon was most thoroughly studied experi- 
mentally in cubic ferromagnetics, such as CdCr 2 Se4, EuO 
and EuS. All these materials, when being deficient of 
chalcogen (oxygen) have n-type conductivity. The deficit 
being strong enough, at low temperatures these materi- 
als have temperature-independent metallic conductivity 
or conductivity with small activation energy (a few meV) ; 
when the temperature T approaches the Curie temper- 
ature T c electroresistivity rises sharply. The rise can 
amount to four orders of magnitude in CdCr2Se4, j^] six 
orders in EuSi-^ and even twelve orders in EuOi-^ 
|^| |l0| . The same behavior of the resistivity can be seen 
in Eu^La^S @, Eu 1 _ 2 .Gd x O jll| and Eux^Gd^Oi^ 
]l2[ . In the paramagnetic region there are two types of 
behavior: in CdCr2Se4 and EuOi_ x (xj.4%)the resistiv- 
ity has large activation energy (~ .2 eV and ~ .3 eV 
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respectively) in the whole temperature range. At the 
same time the resistivity of EuSi-^ and of EuO and EuS 
being doped by a three- valent rare earth metal (and also 
of Eui-zGd-rSe, which at x > .01 is ferromagnetic ||) 
returns at hight temperatures to the initial quasi-metallic 
behavior. 

The common feature of all these materials is a strong 
interaction between electrons and spins of the magnetic 
d- and f-ions. It is commonly accepted now that the elec- 
trical of these materials are determined by the interaction 
of electrons with the magnetic subsystem. Pre-existing 
explanations of the above-mentioned resistivity behavior 
used the following two approaches: 

1. magnetopolaronic approach (in various modifica- 
tions), which was used for the description of metal- 
insulator transition in EuOi-^; 

2. concentrational approach, when it was supposed 
that the resistivity rise is due to the movement of the 
conductivity band edge with the temperature relative to 
the fixed impurity levels. 

We describe the conductivity behavior in all above- 
mentioned materials in a novel approach, which is based 
upon the ideas and methods of localization theory. The 
differences in the resistivity behavior between materials, 
as it will be shown later, are due to the different relations 
of the band width to the spin exchange splitting and also 
to the different electron- impurity interactions. 



III. HAMILTONIAN AND THEORETICAL 
FORMULATION 

Hamiltonian of the electron is 



l.l'a 



(1) 



where V(1T) is the hopping integral, I is the s-d(f) ex- 
change integral, S; is spin localized on / site, cj_, ci„ are 
electron creation and annihilation operators, and V^_i is 
the potential of electron-impurity interaction. (We sup- 
pose that magnetic and crystallographic lattices coincide 
and that there is one electron orbit per site). 

It is convenient to separate out two limiting cases |l| 

1. narrow conduction band W <C 2\I\S, where W is 
the band width; 

2. wide conduction band W ~> 2|/|5. 

As it was shown in band calculations, the conduction 
band in CdCi^Sezi, formed by p-d7 orbitals (3d7 Cr + 
4p Se) is narrow , and p-de conduction bands in EuX 
are wide Jig ], It is worth noting however, that strong 
inequalities for these materials do not hold. 

In Part I of our paper we shall be interested in the nar- 
row band case. In this case it is convenient to consider 
Hamiltonian (pi) in the representation which diagonalizcs 



the main term - s-d exchange. We restrict ourselves to 
the quasiclassical limit 215* 3> 1, when we can consider 
spins as vectors with fixed length, 



Si = Sm 



(2) 



The transformation of the initial creation and annihila- 
tion operators c[ a , ci a into new ones C\ G , Ci„ is equivalent 
to the quantization of the electron spin in the local ref- 
erence frame, the OZ axis being parallel to n. We have 



Ci+ 



Vl 



-viexp{icj)i) uiexp(-i(f>i) J \ c/_ 



cu 



, (3) 



where ui — cos(0;/2),u/ = sin(0//2) exp(— i<f>i). Hamilto- 
nian (|l|) after the transformation (||) takes the form 

H e = ]T v aa \i,i')c} a c lla , + 

< IV > 
aa' 



Y,si aa 'cj a c la , + V e -i, 

laa' 



(4) 



where matrices (in spin space) V(l, V) and / are 



V(l,l') = V 



UlUl 

uivp 

I 

-/ 



- viv v 
vfuv 



— UiVf 



Y VlUv 

■ v\v v 



(5) 



and nearest neighbors approximation is used. 

The conduction band being narrow, we must retain 
in Eq. (^|) the creation and annihilation operators with 
a = sgnl = +1 only. The Hamiltonian then takes the 
form 



E 

<w> 



Vn>Cl + C l/+ + V e -i, 



V ltV = V(n l ,n l >)=V 



CDS I — I COS 



sm|- 



— | expi( 



(6) 



When there is a complete order (T — 0) 0i = 0. When 
T > the angles 8 and 4> and hence the Hamiltonian 
H e are the functions of time. Luckily for us there is a 
small parameter ratio of typical fluctuations frequency to 
electron energy ~ kT c /W. Then the adiabatic approxi- 
mation is applicable, and what we need to solve is static 
problem with the random Hamiltonian H e (t). Hence we 
obtain the disordered system with an off-diagonal dis- 
order |l(|. We suppose also that the influence of the 
electron-impurity interaction on the conduction electrons 
is much less than the influence of magnetic disorder, and 
completely neglect the former. 

The analysis of Hamiltonian (|^) is conveniently started 
with the density of states. 



2 




FIG. 1. The temperature dependence of the density of states in 
doped ferromagnetic semiconductor (the lower local band). Ed 
is the impurity level, £ is the chemical potential, Eq is the band 
edge in the effective medium approximation, ±|-E C | are the mobility 
edges 



IV. THE MOBILITY EDGE CALCULATION 

Let us turn to Hamiltonian (g). As is well known, the 
problem of the mobility edge calculation is equivalent 
to the calculation of the convergence radius for a non- 
averaged Green function self-energy perturbation series 
expansion. The site diagonal representation of matrix 
elements of the Green function being written in the form 



Gu(E) = [E — Mi(e)]~ 



(10) 



we obtain for the self-energy Mi(E) the series 
M t (E) = I £ V lh V hl + ^ V ll2 V hh V hl + ... 



When some approximation of the effective medium is 
used, the influence of the disorder consists in the tem- 
perature dependent narrowing of the band. Such a con- 
ception was used to explain the temperature dependence 
of the electroresistivity of CdCr2Se4 ||. It was supposed 
that the electric current is carried by the electrons acti- 
vated from fixed impurity levels, created by V e -i to the 
red-shifted band edge. For the electrical resistivity the 
following formula was used 



p(T) = po exp 



A_ 

kf 



(7) 



where the activation energy A is defined by 

A = E -E D , (8) 

where Eq is the band edge and Ed is the impurity level. 
There exists a number of facts, however, which do not 
agree with this picture. In CdCr2Se4 above the temper- 
ature T ~ 200 K a large blue shift of the band edge is 
observed Q , but there is no activation energy change in 
this temperature range. Also a similar resistivity transi- 
tion takes place in n-CdCr2S4 |17[| , where the band-edge 
shift is blue at all temperatures g . What is more impor- 
tant, according to modern conceptions for conductivity 
in disordered systems, it is the mobility edge E c which 
is important. We remind that the mobility edge is the 
energy which separates the localized from the extended 
states. For the calculation of the resistivity it is possible 
to use Eq. (§), ||] but Eq. (§) should take the form 



A = E r 



c, 



(9) 



where C is the chemical potential. As in Ref. || we shall 
consider £ as temperature independent. Than our main 
task is to calculate E c . 



Every term in Xl corresponds to a definite path with L 
steps starting from the site l\ and ending at it. To each 
step from a site I to one of its nearest neighbors V there 
corresponds the factor Vw and the passage of every inter- 
mediate site produces multiplication by the factor l/E. 
@ 

It is worth noting that the analysis of the series (^lj) is 
much easier than that of the series involved in the case 
of the traditional Andersen model with diagonal disorder 
|l6| . It is due to the fact that the expectation value of ev- 
ery term in Eq.(|lT|) is finite and, hence, the convergence 
of this series with probability one is equivalent to its con- 
vergence in some norm. Indeed we can use the theorem 
according to which the series of independent random vari- 
ables, say Y]t Xl, , converges with probability one if and 
only if the series J^l < ^ L > an< ^ TIl < -^L > con " 



verge, where the brackets mean the average (in our 
case over the distribution H e (t) , that is due to the ergodic 
hypothesis the average over all spin configurations). Let 
us analyze now the averaged series. We introduce for- 
mally a small parameter l/z, where z is the number of 
nearest neighbors, and then calculate the averages in the 
leading approximation with respect l/z. Since the ratio 
of the number of intersections to the total number of sites 
for a typical path is of the order of l/z p (p > 0), we shall 
on averaging neglect the intersections at all, i.e., when 
calculating < Vi ± i 2 . . . Vi L i x > we consider all the sites in- 
volved as different ones, though the number of paths of of 
the given Length L being approximated by z L Note that 
our analysis deviates from the traditional one for the case 
of the off-diagonal disorder (see e.g. Ref. Neglect- 
ing the correlations due to self-intersections of the paths 
we, however, take into account the correlations created 
by magnetic short-range order, which is more important 
in the problem. 

We are interested in the convergence radius of the se- 
ries 
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<M{E) >=EY^^M L 



(12) 



where the following notation is introduced: 
M L =« V h i 2 V hh . . . Vi Lh >>= 
< J dni . . .dn L V(ni,n 2 ) . . . 
^(n L ,ni)P(ni, . . .n L ) > L 



(13) 



with P(ni,...rii) being the probability density of a 
given spin configuration and the brackets < ... >l de- 
note the average over all paths of a given length. It is 
worth noting that in our approximation the convergence 
of the series J2l < -%l > is equivalent to the convergence 

Let us consider two limiting cases: 
1. T = 

In this case the effective hopping integral does not fluc- 
tuate and the series take the form 



< M{E) >=E^ 



z Lyl 

~1F 



(14) 



and its convergence radius \E C \ = zV naturally coincides 
with the exact band edge (the center of the band corre- 
sponds to E = 0). 
2. T = oo. 

In this case all the paths of a given length are iden- 
tical, all the spin configurations equally participate and 
averaging is easy to make up 



Mr 



1 



dni . . . dn L V(n 1 ,n 2 ) . . . V{n L ,m) (15) 



(4tt) l 

This can be rewritten in the form 
M L = Tr{T/ L }, 



(16) 



where V is an operator with the kernel V(n, n')/47r. 
Then the convergence of the series (|l^) is equivalent to 
the convergence of J2l( z ^) L I E l , where A is the largest 
eigenvalue of the operator V . Thus for the mobility edge 
we obtain 



\EJ = ZX 



(17) 



The eigenvalues of V are easy to obtain due to its factor- 
izability, and the result is 



\E r 



zV 



(18) 



The later case prompts the convergence analysis 
method for arbitrary temperature. Let us use for the spin 
configuration probability the chain superposition approx- 
imation due to Kirkwood |22]| 



P(m,. 



■n L ) = 



P(n 1 ,n 2 )...P(n L ,n 1 ) 
P( ni )...P(n L ) 



(19) 



where P(ni,n') is the pair probability and P(n) is the 
one-site probability. In this approximation all the paths 
are again identical and we again obtain Eq.(p7|), only 
now A is the largest eigenvalue of the operator V with 
the kernel 



V(n,n') = 



V(n,n')P(n,n') 
[^(n)P(n')] 



1/2 



(20) 



Due to the variational principle A is the maximal value 
of the functional 



*"{/(*)} = 



J J r(n)V(n,n')f(n')dndn' 
^/l/(n)| 2 rfn 



(21) 



For the estimation of A we may use the trial function in 
the same form as in the molecular field approximation. 
In this approximation P(n, n') = P(n) • P(n') and the 
operator V is again factorizable. The eigenfunction of 
V corresponding to the larger of the two eigenvalues has 
the form 



/(n) = [P(n)] 1 /2 C os(e/2). 



(22) 



Using this trial function in the variational principle we 
obtain 



IP 



zV 1 + 2 < n z > + < n ■ ni > _ 
2 1+ < n z > 

zV 1 + 2<S Z > [S+ < S ■ Si > /S 2 
~2 1+ < S z > JS ' 



(23) 



Hence the result is the mobility edge expressed in terms 
of magnetization and pair correlation function of the 
nearest- neighbor spins. Our method of mobility edge 
evaluation can be generalized also to the case of arbitrary 
relation between IS and W |2j| . For arbitrary tempera- 
ture the calculation of the mobility edge is equivalent to 
solving the cumbersome integral equation. For T 3> T c , 
however, this equation can be easily solved; the mobility 
edge then is 



\E r 



W 

T 



is 



(is? 
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1/2 



(24) 



V. COMPARISON WITH EXPERIMENT 

Though Eq.(|2^) is obtained for classical spins we'll use 
it to explain the temperature dependence of the resistiv- 
ity of CdCr2Se4 (Fig. 2). The temperature dependence 
of the resistivity at low temperatures (T < 50 K) shows 
that the chemical potential lies in the close vicinity of 
the band edge. At those temperatures the resistivity is 
defined by the exact position of £ relative to the band 
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edge; at higher temperatures when E c moves upward, 
the dependence of p on T is determined mainly by the 
shift of E c and for the chemical potential we may use 
the approximation £ = — zV. Inserting into Eq. (R) the 
expression for the mobility edge we obtain 



Po exp 



- < Sp ■ Si > /S 2 t W 
1+<S Z > /S \AkT 



where W 
for the correlator < So 



(25) 



Using 



2zV is the conductivity band width. 

Si > the results of Ref. |24| , and 
for < S z > the molecular field approximation (S = 3/2), 
and choosing the appropriate values of W and po, we 
obtain the plot presented on Fig. 3. 
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FIG. 2. Electrical resistance of Cdo.9gIno.oiCr2SE4 (reproduced 
from Ref. . 

FIG. 3. Electrical resistance calculated on the basis of Eq. (||). 



The chosen value for the band width is W = 0.8 eV, 
which well agrees with other estimations jn§[l4| ■ We can 
deduce that Eq. (^5|) can claim for a quantitative de- 
scription of the experimental dependence. This formula 
also explains the colossal negative magnetoresistance of 
the materials under consideration, as the magnetic field 
strongly influences the spin correlators in the exponent, 
especially near T c . 

In Refs. PJ25| the data for high resistivity samples of 
CdCr2Se4 was presented. In this case the resistivity as 
the function of temperature has no maximum, but the 
dependence of logp on 1/T changes its slope after the 
transition to the paramagnetic region. Such a behavior 
is also described by Eqs. (|?]) and (^3|) with the chemical 
potential ( < —zV. 



VI. DISCUSSION (10 YEARS AFTER) 

It is expedient to discuss the applicability of the ap- 
proach of this ancestor paper to new classes of materials 
and its place among more recent theoretical approaches. 



Our principal point was that strongly interacting sys- 
tem of electrons and magnetic ions (s-d exchange model) 
in classical approximation for the ions spins and narrow 
electron band, reduces to a random-bond model with off- 
diagonal disorder depending on an instant spin configu- 
ration. This double-exchange (DE) Hamiltonian (Eq. (||) 
in the present manuscript) has been reopened quite re- 
cently 

The model, we believe, can be applied with min- 
imal changes to manganite pyrochlores discovered re- 
cently p7pS| ]. Its application to high- carrier concen- 
tration CMR compounds such as manganite perovskites 
demands discussions and probably additional ideas. For- 
mally, in this approximation the difference between, say, 
LaSrMnOs and CdCr2Se4 (except for carriers concentra- 
tion) is that in the latter the doping electrons appear (due 
to Se vacancies) in the empty d-band, whereas in the for- 
mer the doping holes appear in the filled d-band (due to 
capture of d-electrons by Sr atoms). According to our 
concept, it is the movement of mobility edge, with tem- 
perature and magnetic field, relative to the Fermi level 
that explains the resistivity peak and CMR effect. It is 
concordant with the recent conclusion of Varma |^9| who 
suggested the decisive role of the above spin-dependent 
disorder in metal-semiconductor transition and CMR ef- 
fect in manganite perovskites. 

The viewpoint that spin disorder alone can not local- 
ize carriers at temperatures close to T c and in paramag- 
netic state has been argued by Millis et al. |30| on the 
basis of a perturbational calculation. Later this group 
presented very subtle arguments concerning probability 
distribution of spin-dependent hopping integral |nj in fa- 
vor of their conclusion. The matter is that, as with every 
random-bond model, DE Hamiltonian keeps a large por- 
tion of the states at the band center delocalized irrespec- 
tive of how the bond fluctuations are strong. Therefore, 
it may turn out that the weight of the localized states 
be insufficient to embed large number of holes (0.2 to 
0.3 per cell). After numerical simulation with the use of 
one-parameter scaling localization theory |32j it has been 
found that this is really the case for pure DE mechanism. 
Note that for SC lattice our result of 1988 E c = -W (see 
manuscript above) is in a reasonable agreement with the 
new more exact E c = — 3.55U p2^ . 

Thus, an agreement with experiment for manganites 
may be achieved by introducing some additional mech- 
anism. Millis et al. |3l[] developed and advocated Jahn- 
Teller polaron mechanism to explain the resistivity peak 
and CMR. On the other hand, in 1988 we proposed the 
mechanism based on the pinning of the Fermi energy by 
the donor levels, due to very weak overlap of the vacancy 
states with magnetic d-states in CdC^Se^see present 
paper). Recently our mechanism has been recognized 
and applied to LaBaMn0 3 by Bebenin et al. Q. 
In the case of manganite perovskites one may speculate 
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that the Fermi level is more or less fixed by impurity 
band formed by Sr, Ba, etc. acceptor states. The band 
emerges in Mn d-band at low temperatures and holes 
are mobile because the d-band is delocalized. With in- 
creasing the temperature the mobility edge shifts towards 
higher energies and holes remain in localized part of the 
band. 

Recently an approximate pinning of the Fermi level 
has been found in a model where diagonal disorder is 
added to DE hopping Hamiltonian [34] ]. Numerical sim- 
ulation using one-parameter scaling localization theory 
have shown that the states at the Fermi energy are lo- 
calized in PM region if the intensity of diagonal disorder 
is as intense as about the bandwidth fjij . Though this 
model is very attractive, the physical origin of such strong 
fluctuation of Mn d-levels is not evident. Sheng et al. |3J] 
attribute it to potential fluctuations due to random sub- 
stitution of La by Sr, Ca etc. using experimental data 
on resistivity of polycrystalline samples. However, con- 
trary to the assumption of |34| there exist single crystals 
with zero-temperature resistivity smaller than 10 -4 f2-cm 
which do exhibit large resistivity peak [f35f . This is an 
evidence in favor of weak substitutional disorder, and 
impurity band formation. 

The mechanism of conduction in PM region is not es- 
tablished until now but it lies out of scope of this discus- 
sion. We hope to return to discussion on this problem in 
a forthcoming paper. 
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